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ABSTRACT 

Observations of quasi-periodic oscillations (QPOs) in X-ray binaries hold a key to understanding 
many aspects of these enigmatic systems. Complex appearance of the Fourier phase lags related 
to QPOs is one of the most puzzling observational effects in accreting black holes. In this Letter 
we show that QPO properties, including phase lags, can be explained in a framework of a simple 
scenario, where the oscillating media provides a feedback on the emerging spectrum. We demonstrate 
that the QPO waveform is presented by the product of a perturbation and a time delayed response 
factors, where the response is energy dependent. The essential property of this effect is its non-linear 
and multiplicative nature. Our multiplicative reverberation model successfully describes the QPO 
components in energy dependent power spectra as well as the appearance of the phase lags between 
signal in different energy bands. We apply our model to QPOs observed by Rossi X-ray Timing 
Explorer in BH candidate XTE J 1550-564. We briefly discuss the implications of the observed energy 
dependence of the QPO reverberation times and amplitudes to the nature of the power law spectral 
component and its variability. 

Subject headings: accretion, accretion disks — black hole physics — X-rays: binaries — stars: individual 
(XTE J1550-564) 


1. INTRODUCTION 

Galactic X-ray Black Hole (BH) candidates exhibit 
spectacular evolution of timing and spectral behavior 
(Remillard & McClintock 2006; Belloni 2005). Strong 
variability is usually observed in states, dominated by 
non-thermal emission (i.e. hard and intermediate states). 
The variability properties strongly correlate with over- 
all spectral state evolution (Klein- Wolt & van der Klis 
2008). A common property of the variability is nearly 
periodic flux modulations with frequency between 0.01 
and 20 Hz generally referred to as low frequency quasi- 
periodic oscillations (QPOs). The exact origin of QPOs 
in BH systems is still under debate. In this Paper we 
present a new important evidence which allows us to con- 
strain the QPO models and to advance in understanding 
of this phenomenon and as well as the overall BH phe- 
nomelogy. 

QPOs are observed as one or more relatively narrow 
peaks in Fourier power density spectra (PDS). QPO sig- 
nal is shown to be strongly coupled with the power law 
part of the energy spectrum (Shaposhnikov et al. 2011; 
Sobolewska & Zycki 2006). Probably the most intrigu- 
ing aspect of the QPO phenomenology is the behavior of 
Fourier phase lags, i.e. differences in phases calculated in 
two different energy ranges. High time resolution data on 
brightest BH transients, provided by Rossi X-ray Timing 
Explorer ( RXTE) allowed to study this effect in detail. 
Phase lags observed during the initial intermediate state 
of 1998 event observed in XTE J1550-564 have shown 
complex behavior (Cui et al. 2000; Remillard et al. 2002). 
Particularly remarkable is the fact that the fundamental 
QPO and its first harmonic exhibit lags of opposite signs. 
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While the fundamental has shown the negative (or soft) 
lags, the lag in the first harmonics was hard for most ob- 
servations (conventionally we define positive or hard lags 
when signal is delayed at harder energies and vise versa 
for negative or soft lags) . Similar behavior was observed 
during 1999 outburst from XTE J1859+226 (Casella et 
al. 2004) and in persistent BH candidate GRS 1915+105 
(Lin et al. 2000). The conclusion of these studies was 
that either the lags behavior is due to different physi- 
cal origin of QPO harmonics or some unknown physical 
mechanism lead to a particular QPO waveform depen- 
dence on energy, which is well reproduced in many obser- 
vations. However, no QPO model have met the challenge 
to account the observed phase lag behavior. 

In this Letter we present new and extremely impor- 
tant paradigm which allows much better understand- 
ing of QPO behavior and put a stronger constrains on 
its physical origin. First, we present a QPO wave- 
form paxametrization which provide proper description 
of QPO properties both in term of Fourier amplitudes 
and phases. In fact, we develop the first analytical model 
for QPO phase lags which allows direct modeling of the 
phase lags seen in observations. Secondly, based on the 
results of our model, we draw conclusions on the physi- 
cal nature of the observed oscillations and the oscillating 
media. We tentatively call our model the ”QPO rever- 
beration model” . We do not introduce any new physical 
paradigm. Instead, acting on the premises of the stan- 
dard truncated disk/hot inner corona model, we identify 
a specific way, in which all major observed QPO prop- 
erties are naturally produced. The main idea behind 
the ”QPO reverberation” model is that the inner part 
of the accretion flow acts as a non-linear forced oscilla- 
tor where both perturbation and response signals enter 
the observed signal in multiplicative manner. While the 
perturbation is hydrodynamical in nature, the response 
factor reflects the spectral shape reaction to the changes 
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in oscillating media properties introduced by the initial 
perturbation. 

The QPO reverberation model matches the data very 
well, allowing to convert a complicated form of the phase 
lags into well understood physical delay time. The re- 
sponse signal is always delayed with respect to the per- 
turbation, as expected in case of forced oscillations. The 
main conclusion which we draw based on our analysis 
is that multiplicative reverberation is a primary mech- 
anism leading to the non-sinusoidal nature of the QPO 
waveform and defining the behavior of QPO phase lags. 

In the next section we describe the model and our 
RXTE data analysis. In §3 we briefly discuss the impli- 
cations of our results to different QPO scenarios. Con- 
clusions follow in the §4. 

2. REVERBERATION EFFECT IN QPO 
2.1. Model 

We parametrize the QPO signal by the following wave- 
form: 

s(t) oc (l + ae'^cosa/o^Cl + ^^cos^oCt-td))). (1) 

We assume that an oscillation process in the system 
signified by QPO leads to a time delayed response of the 
spectral form of the outgoing radiation. The motiva- 
tion for this functional form of the signal is most easily 
understood In terms of the power law shape of the av- 
erage energy spectrum NE~ a . If we assume that the 
normalization N and the index a experience exponen- 
tially decaying oscillations, where oscillations in index 
are delayed with respect to the oscillations in the nor- 
malization then, after expanding the power law in Tay- 
lor series keeping only the first term, we arrive to the 
expression (1). The Fourier transform of the waveform 
and its products (i.e. power spectrum, phase lags) are 
presented in Appendix. 

2.2. Application to data 

We apply our model to the data collected with RXTE 
during the brightest BH transient event observed from 
XTE J 1550-564 (Homan et al. 2001; Cui et al. 2000; Sha- 
poshnikov & Titarchuk 2009). We use a combination of 
Binned and Event data modes to extract 7.8125 millisec- 
ond time resolution lightcurves in channels covering 1.5- 
3.5, 3.5-5.0,5.0-6.0, 6.0-7.5, 7.5-12.0, 12.0-14.6, 14.6-19.3 
and 19.3 - 50.0 keV energy ranges. We compute Fourier 
transforms for consecutive intervals of 64.0 sec for each 
lightcurve. We then compute corresponding power spec- 
tra and phase lags as a complex argument of the cross 
spectrum using the first channel range as a reference. We 
rebin these data products using 1.02 logarithmic rebin- 
ning factor. In Figures 1 and 2 we show representative 
power spectra and phase lags for two observations, 30188- 
06-05-00 (Observation 1) and 30191-01-30-00 (Observa- 
tion 2). The fundamental QPO during these observations 
appeared at 2.35 Hz and at 6.5 Hz respectively. The first 
observation is taken during a harder spectral state, i.e. 
earlier in the outburst. It showed positive lags for both 
QPOs, while the second observation shows negative lag 
at the fundamental. For each observation we combine 
eight PDSs and seven corresponding phase lags spectra 
in a joint model fit. Using the QPO waveform described 
by Eq. 1 we model QPOs in the power spectra and cor- 
responding phase lags between signals. The parameters 



FlG. 1. — RXTE data from observation of XTE J1550-564 (Ob- 
sID 30188-06-05-00). The bottom panel shows phase lags between 
signals in soft and hard energy ranges. The middle panels show 
PDSs for the two energy ranges. The data are fitted with the re- 
verberation model, shown as solid color lines (red for the soft range 
and blue for the hard). For the PDS fits the dashed lines show con- 
tribution of Poisson and broad-band noise components. The top 
panel shows the QPO waveforms for signal for both energies the 
soft (green) and the hard energy ranges as given by the model fit 
parameters. 

describing perturbation signal a and A are the same for 
all waveforms, while the parameters describing spectrum 
reverberation, i.e. its amplitude b and time delay td (or 
phase delay <j> — 27tcjo td) are individual for each energy 
range. The data and the model fits for two energy ranges 
are presented in Fig. 1 and 2. 

Some systematic is seen in the phase lags fits in the 
vicinity of the first harmonic of the QPO in power spec- 
tra. It would be naive to expect that the parametriz ac- 
tion using the expression (1) will describe every aspect 
of the data. Rather than the exact QPO waveform, the 
presented model provides a basic framework for develop- 
ment of the future models which capture such effects as 
contribution of aperiodic PDS component in the phase 
lags, effects of statistical noise, etc. Taking this fact, 
the model fits the data extremely well capturing major 
behavior of the phase lags and the exact form of the 
fundamental QPO, which is energy dependent and assy- 
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Fig. 2. — The same ss Fig 1 for RXTE observation 30191-01-30- 
00 which is performed later during the outburst. Note that for this 
case, as opposed to the case presented in Fig 1 the phase lags at 
the fundamental and the first harmonic have different signs. 

metrical. In Figures 3 and 4 we show the dependence of 
the reverberation times td and amplitudes b on energy. 
Below we discuss the implications of these results for the 
physical picture of QPO production in BH candidates. 

3. DISCUSSION 

The results of the presented model fits to the data show 
that the reverberation scenario properly describes the 
QPO properties. Apparent non-linear multiplicative na- 
ture of the effect and its energy dependence naturally ex- 
plains a number of observational effects including shape 
of the QPOs in PDS, energy dependence of the QPO 
powers, evolution and appearance of the phase lags. 

Probably even more important result of our analysis 
is the fact that negative Fourier time lags observed in 
many cases for the fundamental QPO are of a purely 
mathematical nature and should not be treated as physi- 
cal times. The reverberation model properly convolves a 
complicated behavior of Fourier phases into a one, well 
understood, physical time td which can then be consid- 
ered in terms of physical models. 

Another important result is the nature of the observed 
reverberation and energy dependence of the reverbera- 



Fig. 3. — Energy dependence of the reverberation time td for the 
Observation 1 (blue points, a = 1.30 ± 0.06, A — 0.56 ± 0.01) and 
Observation 2 (red points, a = 1.0 ± 0.1, A = 2.07 ± 0.03). The 
logarithm f(E) — Alog(^) 4- B fit is shown by black solid line 
for the first observation and dashed line for the second one. The 
time delays for the first observations clearly follow the logarithm 
dependence while the data for the Observation 2 is consistent with 
the constant delays. 

tion time and amplitude . It is clear that the perturba- 
tion actively changes the condition in the system. Geo- 
metrical QPO models which are based on a mechanical 
rotation or precession (Stella & Vietri 1998; Schnittman 
et al. 2006, see e.g.) lack this active ingredient. Rever- 
beration effect strongly favors QPO models which involve 
oscillations affecting the state of matter which then is re- 
flected in the form of the emerging spectrum. This group 
of models includes transition layer (corona) oscillations 
(Titarchuk et al. 1998) and gravitationally trapped disk 
oscillation modes (Wagoner et al. 2001). Logarithmic de- 
pendence of the time delays for the hard state (see Fig. 
3 for Observation 1) strongly favors accretion induced 
oscillations of the Compton corona as model for QPOs. 

In the truncated disk scenario a thin accretion disk 
extends down to some inner radius where it superseded 
by a optically-thin geometrically thick configuration or 
corona (see e.g. Tomsick et al. 2009). In this scenario, 
the thin disk is responsible for the soft thermal part of 
the spectrum, while the corona produces a non-thermal 
power law part of the energy spectrum. Variation in the 
inner disk radius, caused by surges of accreting matter, 
introduce perturbation in the corona and modulate its 
physical parameters, i.e. its density, temperature, optical 
depth, etc. This leads to temporal variations of the spec- 
tral shape of the escaping radiation. The time delay be- 
tween the perturbation and the spectral response is set by 
sound speed in the corona. The logarithmic dependence 
of time delays with energy observed for Observation 1 
(hard state) indicates (see Figure 3), according to ex- 
pectations, that Comptonization play significant part in 
spectral formation during these states. Namely, multiple 
upscattering events experienced by photons cause wave- 
like structures to propagate from the low energies toward 
the higher part of the time dependent spectrum. On the 
other hand, during the soft intermediate state (Obser- 
vation 2), the reverberation time is flat with respect to 
energy, indicating that the spectrum reacting to the per- 
turbation at the same time at all energies. Moreover, the 
reverberation amplitude for this observation is well rep- 
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Fig. 4. — Energy dependence of the reverberation amplitude b for 
the Observation 1 (blue) and Observation 2 (red). The Observation 
2 data is fitted with f(E) = A log(l /E)-\-B which is shown by black 
solid line. 

resented by b oc log(l/jE), in agreement with variability 
being produced by a pure pivoting power law. This can 
not be facilitated through the thermal Comptonization 
and will require either synchrotron or non-thermal (such 
as bulk motion, see e.g. Titaxchuk&Shaposhnikov 2010) 
Comptonization mechanisms to be invoked. This subject 
is beyond this Letter and will be addressed in more detail 
in upcoming publications. 

Maccarone et al. (2011) studied QPOs observed in 
GR3 1915+105 by means of bispectrum and found evi- 
dence for strong non-linear coupling between QPO har- 
monics and the broadband noise and considered model 
involving non-linear oscillator. Our findings also strongly 
suggest that the inner accretion flow region acts as a non- 
linear forced oscillator, which, by means of combination 


of radiation mechanisms such as thermal/non-thermal 
Comptonization and synchrotron radiation, modulate 
the emerging spectrum of the system, leading to a non- 
linear nature of the QPO signal. The waveform presented 
by expression (1) predicts that on time scales shorter 
than the exponential decay time 1/A both the average 
flux and variability are proportional to e~ Xt and, there- 
fore, are linearly related to each other. Rms-flux relation, 
well established in broad variability from BH sources, was 
recently found in QPOs from XTE J1550-564 by Heil et 
al. (2011). A specific feature of the relation is that it 
switches from positive to negative above 5 Hz. Heil et 
al. (2011) show that scaling of the QPO power to the 
power of underlying broad band noise corresponding to 
the QPO frequency reinstates the linear rms-flux rela- 
tion, in excellent agreement with the expectation that 
the noise variability is associated with the driving force 
for non-linear oscillator associated with QPO signal. 

4. CONCLUSIONS 

We present a new observational effect in QPOs. We 
identify the time delayed reverberation in QPO signal as 
mechanism facilitating the phase lag behavior observed 
in bright galactic BH candidates. The QPO reverbera- 
tion model successfully describes the QPOs in energy de- 
pendent power spectra and the corresponding phase lags. 
The presented results strongly point to oscillation of the 
Comptonizing and, most probably magnetized, corona as 
the origin of the QPO signal from accreting BHs. The 
QPO reverberation parametrization presents, for the first 
time, a consistent description of the QPO signal both in 
Fourier amplitude and phase domains, providing an ex- 
citing avenue to finally resolving the QPO nature. 


APPENDIX 
QPO WAVEFORM 

Let us consider the waveform described by the expression: 

s(£, E) oc (1 + ae~ Xt coswo£)(l + b(E)e~ Xt cos(cjo* — <£(#))), 


where b(E) and <f>{E) are functions of energy E. 
The Fourier transform of s(t) is given by 


S&,E) = - 
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(A2) 


The QPO power spectrum and phase lags between signals at energies E± and E 2 are calculated in straightforward 
manner, i.e.: 

P(E) = \S(E )\ 2 (A3) 


and 

correspondingly. 


SP(E U E 2 ) = Arg (SfaEh) * S(u,E 2 )) 


(A4) 


REFERENCES 


Belloni, T. 2005, Interacting Binaries: Accretion, Evolution, and 
Outcomes, 797, 197 

Casella, P., Belloni, T., Homan, J., Sc Stella, L. 2004, A&A, 426, 
587 

Casella, P., Belloni, T., & Stella, L. 2005, ApJ, 629, 403 
Cui, W., Zhang, S. N., Sc Chen, W. 2000, ApJ, 531, L45 
Heil, L. M., Vaughan, S,, Sc Uttley, P. 2011, MNRAS, 411, L66 


Homan, J., Wijnands, R., van der Klis, M., et at 2001, ApJS, 
132 377 

Klein-Wolt, M., & van der Klis, M. 2008, ApJ, 675, 1407 
Lin, D., Smith, L A., Liang, E. P., Bottcher, M. 2000, ApJ, 543, 
L141 

Maccarone, T. J., Uttley, P., van der Klis, M., Wijnands, 

R. A. D., Sc Coppi, P. S. 2011, MNRAS, 413, 1819 


5 


McClintock, J. E., Remillard, R. A., Rupen, M. P., et al. 2009, 
ApJ, 698, 1398 

Remillard, R. A., Sobczak, G. J., Muno, M. P., &; McClintock, 

J. E. 2002, ApJ, 564, 962 

Remillard, R. A., & McClintock, J. E. 2006, ARA&A, 44, 49 
Shaposhnikov, N., & Titarchuk, L. 2009, ApJ, 699, 453 
Shaposhnikov, N., Swank, J. H., Markwardt, C., & Krimm, H. 
2011, arXiv: 1103. 0531 

Schnittman, J. D., Homan, J., & Miller, J, M. 2006, ApJ, 642, 420 
Sobolewska, M. A., Zycki, R T. 2006, MNRAS, 370, 405 
Stella, L., & Vietri, M. 1998, ApJ, 492, L59 
Sunyaev, R.A. & Titarchuk, L.G. 1980, A&A, 86, 121 


Tomsick, J. A., Yamaoka, K., Corbel, S., et al. 2009, ApJ, 707, 
L87 

Tagger, M., & Pellat, R. 1999, A&A, 349, 1003 
Titarchuk, L., Lapidus, I., & Muslimov, A. 1998, ApJ, 499, 315 
Titarchuk, L., &: Shaposhnikov, N. 2010, ApJ, 724, 1147 
Vignarca, F., Migliari, S., Belloni, T., Psaltis, D., & van der Klis, 
M. 2003, A&A, 397, 729 

Wagoner, R. V., Silbergleit, A. S., & Ortega^ Rodriguez, M. 2001, 
ApJ, 559, L25 



